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An  expression  is  found  for  the  flux  due  to  earthshine 
incident  at  any  point  on  a  spherical  target  as  a  function  of 
target  height  above  ground.  Curves  showing  the  variation  of 
reflected  signature  (or  view  factor)  with  observer  position  and 
target  height  are  presented  for  the  case  of  a  diffuse  target 
surface.  An  expression  is  given  for  the  view  factor  in  the 
specular  case,  which  is  shown  to  be  independent  of  observer 
position. 
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I. 


INTRODUCTION 


An  important  component  of  the  infra-red  (IR)  signature  of 
a  target  outside  the  atmosphere  is  the  earthshine  it  reflects. 
This  report  will  present  curves  showing  the  amount  of  earthshine 
reflected  by  an  idealized  spherical  target  towards  an  IR  sensor 
as  a  function  of  target  altitude  and  sensor  position  for  the  two 
limiting  cases  of  diffuse  and  specular  target  surfaces.  A 
cloud  of  convex  particles  of  arbitrary  shape  and  random  orienta¬ 
tion  will  exhibit  the  same  behavior  as  a  sphere  of  the  same 
surface  characteristics. 

The  amount  of  energy  incident  on  an  area  element  dA  of  the 
target  in  a  wavelength  band  AX  is 


E  ( AX ) 


/ 

all  dS 
visible 
to  dA 


Ne(AX)  cos0  cos 4’ 


dS 


(1) 


where  E(AX)  is  the  energy  incident  on  dA  (in  Watts) 

Nf(AX)  is  the  radiance  of  the  earth  at  dS  in  Watts  per 
square  meter  per  projected  steradian-assumed 
constant  over  the  earth's  surface. 

0  and  il>  are  the  angles  between  the  line  joining  dA  and  dS 
and  the  local  normals. 

£  is  the  distance  between  dA  and  dS. 


(See  Fig.  1) . 

If  the  bidirectional  reflectance  of  the  material  is  known. 
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s 


then  the  power  reflected  per  steradian  in  the  direction  of  the 
sensor  can  be  calculated  as 


J(AA)  =  f  f 

all  dA  all  dS 
visible  visible 
to  sensor  to  dA 


N_ (AA) cosO 

£j 


cost \>  cos£ 


(12} 

i2 


dS  dA 


(2) 


where  <j>^  and  4>2  are  the  azimuthal  angles  of  dS  and  the  sensor 
measured  at  dA.  To  simplify  calculation,  in  this  report  we 
characterize  the  fraction  of  the  input  power  which  is  reflected 
with  a  single  number  p  {which  is  therefore  assumed  independent  of 
incidence  angle  and  wavelength  over  the  range  of  measurement) . 
Where  this  reflected  power  is  directed  will  depend  on  our 
assumptions  about  the  surface.  A  smooth  surface  will  reflect 
specularly,  i.e.,  the  output  distribution  is  a  narrow  peak 
centered  on  A  rough  surface  will  diffuse  the  output  power 

isotropically  into  the  hemisphere  above  it. 

This  simplification  permits  us  to  split  up  the  integral  into 

two  parts,  and  first  find  the  incident  energy  per  unit  area  at 
any  point  on  the  target  sphere  as  a  function  of  postion  on  the 

sphere,  and  the  sphere's  height  above  ground.  From  this,  the 
total  amount  of  energy  incident  upon  the  target,  and  the  fraction 
of  it  which  is  reflected  toward  a  sensor  in  a  given  direction 
can  be  calculated. 
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II.  DISTRIBUTION  OF  INCIDENT  FLUX  ON  THE  SPHERE 

Tc  find  out  how  much  energy  is  falling  on  any  small  patch 
n  '  the  target  sphere  with  area  dA,  it  is  necessary  to  integrate 
the  radiation  emitted  in  the  direction  of  the  patch  over  all  the 
earth  visible  to  that  part  of  the  sphere. 


Energy  density  at  the  position  of  dA  =  E 


cosijj  cos9 
- 2 —  dS 

visible 

surface  ( 3 ) 

of  earth 


•/here  i  is  the  distance  from  the  sphere  to  dS.  (See  Fig.  2)  . 

This  is  a  very  difficult  problem  for  all  but  the  simplest 
cases,  but  for  an  isotropically  emitting  earth  it  is  permissible 
to  replace  the  earth  by  a  flat  circular  disc  of  the  same  emittance 
subtending  the  same  solid  angle  at  the  sphere.  The  incident  flux 
at  dA  is  then 


// 


cos0  cos  ip  tan0  d0d<f> 


visible 

earth-disc 


where  the  portion  of  the  earth-disc  visible  to  dA  is  described 

by  0  .  <  9  <  3,  and  -p  <  <j>  <  p  as  shown  in  Fig.  3,  where  3 
min  —  —  —  — 

is  the  horizon  angle  and  p  depends  on  x,  the  position  of  dA  on 
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»*•*,*#*_ 


Fig.  2.  Earthshine  incident  on  spherical  target. 
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Three  regions  of  the  sphere  may  be  distinguished,  each  with  a 
different  form  of  expression  for  E.  The  arrangement  of  the 
regions  is  shown  in  Fig.  4(a),  and  their  sizes  as  a  function  of 
height  in  Fig.  4(b). 

Region  1  Contains  those  areas  of  the  sphere  which  can  "see"  all 
the  earth-disc  out  to  the  horizon,  and  are  not  shaded  by  other 
parts  of  the  sphere.  The  range  of  values  for  x  is 

0  <  x  <  j  -  6  (6) 

and  for  this  region 

0  .  =0;  p  =  tt 

mm  r 

since  the  integration  is  over  the  whole  disc.  This  yields 

2 

E  =  NgiT  cos  x  sin  3  (7) 

Region  2  Contains  those  areas  of  the  sphere  which  can  "see" 
some  of  the  earth-disc,  but  are  shaded  by  other  parts  of  the  sphere 
from  seeing  the  remainder.  If  a  part  of  the  earth  can  be  seen 
by  a  given  patch  on  the  sphere,  it  will  have 

cos  =  cos  <p  sin  8  sin  x  +  cos  0  cos  x  >_  0  (8) 

Hence  for  the  range 

TT  TT 

2  -  e  <  x  <  J  +  8 
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1113683  -N 


we  have 


®min  ^an  ^eosi})  tan  xj  '  ^  cos  |tanB  tan 

which  after  considerable  manipulation  yields 


(9) 


E  =  N„{p  cos  x  sin^B-sin  p  sin  x  sin  B  cos  B-tan  ‘L(tan  p  cos  x)  } 


-1 


=  Nt 


cos 


-1 


•1/ _ ■!--  \ 

ItanB  tan  x/ 


2  (22 
cos  x  sin  B-cos  B  \ sin  x-cos  B 


(10) 


(11) 


+  cos  (cos  B/sin  x) 

Region  3  Includes  those  areas  of  the  sphere  which  are  completely 
screened  from  any  view  of  the  earth.  This  is  true  for  the  areas 
which  have 


x  >  2  +  3 


Not  surprisingly,  this  gives 


E  =  0 


(12) 


(13) 


If  the  height  of  the  sphere  above  ground  is  h,  and  the  earth's 
radius  is  R,  then 


sin  s  ■  ra  =  ifer  if  “  ■  h/R 


(14) 


The  function  E'  =  E/N„  is  plotted  in  Fig.  5  as  function  of 
x(in  degrees)  with  h  as  a  parameter.  To  facilitate  comparison 
of  the  angular  dependence  of  E'  at  different  heights,  Fig.  6 
shows  the  same  curves  but  normalized  to  unity  at  x  =  0  (i.e., 


10 


INCIDENT  FLUX  E'*E/N 


0  20  40  60  80  100  120  140  160  180 


POSITION  ANGLE  ON  THE  SPHERE  (d*g) 


Fig.  5.  Incident  flux  as  a  function  of  position  on  the 
sphere  for  various  target  heights. 
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the  curves  are  of  E' (x, 8)/E' (0,6)).  This  format  shows  the  differ 
ent  behavior  in  the  three  regions;  i.e.,  following  the  cosine 

curve  in  region  1,  intermediate  behavior  in  region  2  and  identi¬ 
cally  zero  in  region  3. 

Special  Cases  When  h  is  either  very  small  or  very  large,  the 
expression  for  E  can  be  simplified. 

As  h  0,  8  -*•  tt/2,  and  the  region  2  solution  is  valid  for 
most  of  the  sphere.  Looking  at  this  solution  as  M  tt/2  and 
therefore  p  ir/2  we  find 

E(x,  8  +  tt/2)  =  Neit(1+cos  x)/2  (25) 

At  the  other  extreme,  when  the  earth  is  very  far  from  the 

sphere,  it  acts  like  a  point  source,  and  the  region  2  disappears. 
Hence 

E(x,8  +  0)  =  Ne  ttcos  x  sin28  x  <  ir/2  (16) 

=  0  x  >  tt/2 


III.  TOTAL  POWER 


From  a  knowledge  of  the  incident  flux  on  the  sphere  as  a 
function  of  position,  the  total  power  incident  on  the  sphere  can 
be  found  at  any  height  by  integrating  E  over  its  surface.  In 
practice,  the  form  of  equation  (11)  makes  the  integration  rather 
arduous ,  but  it  eventually  yields  the  result 


P  =  2N£  irr  it  (1-cos  8)  (] 

=  2N_  A  it  (1-cos  8) 

E  p 

2 

where  r  =  the  sphere's  radius  and  Ap  =  Trr  =  the  projected  area 
of  the  sphere. 


IV.  OBSERVED  BRIGHTNESS 


The  infra-red  brightness  of  the  sphere,  as  seen  by  a  sensor 
in  the  far  field,  is  obtained  by  integrating  the  power  reflected 
by  the  sphere  in  the  direction  of  the  sensor,  over  the  hemisphere 
visible  to  the  sensor. 

Case  (a) .  Specular  Surface 

If  the  surface  of  the  sphere  reflects  specularly,  then  its 

spherical  shape  will  diverge  an  incoming  beam  of  radiation  so 

as  to  put  equal  power  in  all  directions,  as  will  now  be  shown. 

Let  a  plane  wave  be  incident  on  a  specularly  reflecting 

sphere  of  radius  r.  (See  Fig.  7) .  To  calculate  the  power 

scattered  by  the  sphere  as  a  function  of  angle  measured  from  the 

direction  of  origin  of  the  plane  wave,  we  find  the  power  incident 

on  an  imaginary  enclosing  sphere  of  radius  R>>r  as  a  function  of 

position.  If  r>>A  the  wavelength  of  the  radiation,  each  (02,  tp p) 

point  on  the  outer  sphere  will  receive  radiation  from  only  one 

2 

point  (91,<});l)  on  the  inner  sphere.  The  small  area  dA^=r  sin0d0d<j> 

on  the  inner  sphere  has  a  total  incident  power  of  P  dA^  cos  0 

2 

Joules/sec  if  P  Joules/m  sec  is  the  power  density  in  the  plane 

wave.  A  fraction  (p)  of  this  is  reflected.  The  corresponding 

2 

area  on  the  outer  sphere  is  dA2=R  sin02d02d<f>2 '  an<^  ky  the  ^aw  of 

specular  reflection  02=20^;  d02  =  2d0^;  <j)2=<J>^.  This  means 
2  2 

dA2=4(R  /r  JdA-^  cos  0.  All  the  power  reflected  by  dA^  falls  on 


dA2,  which  therefore  has  a  received  power  density  equal  to: 


pPdA^  cos  9^  pPr 


dA, 


4R 


independent 

of 

position. 


Hence  a  plane  wave  incident  on  a  specularly  reflecting  sphere 
is  reflected  in  all  directions  equally.  The  radiant  energy 
from  the  earth  can  be  considered  as  made  up  of  an  angular  spectrum 
of  plane  waves.  Each  plane  wave  will  be  reflected  isotropically 
and  so  the  total  brightness  also  will  be  isotropic. 

The  received  power  per  steradian  of  sensor  is  therefore  just 
(the  total  power  reflected  from  the  sphere)  /4tt,  which  is 

J  =  2N.-A  tt  (1-cosB)  p/4ir  (18) 

t  P 

J  =  N_A  p ( 1-cos 8) /2  (19) 

£<  P 

A  va_ue  frequently  used  to  characterize  a  particular  geometry 
is  the  view  factor  F,  defined  by 


F  = 


VAP 


(20) 


For  the  specular  case  therefore, 


F  =  (l-cos8)/2 


(21) 


and  is  independent  of  sensor  position. 


1.1; 
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Case  (b) .  Diffuse  Surface 


If  the  sphere  is  a  diffuse  reflector,  then  the  observed 
brightness  will  depend  on  the  sensor  angle  y  (see  Fig.  8).  The 
radiant  intensity  seen  by  the  sensor  .'.3  obtained  by  integrating 
(E  cos£)  /ir  (since  the  incident  earthshine  is  scattered  into  it 

projected  steradians  )  over  the  half  of  the  sphere  visible  to 
the  sensor  (0£|£|<J).  In  the  previous  section  expressions  for 
E  were  obtained,  but  as  a  function  of  x,  which  is  not  equal  to 
£  unless  y=0  (sensor  below  target).  Expressing  E  as  a  function 
of  £  and  y  in  the  general  case  leads  to  a  form  too  complicated 
fcr  analytical  integration.  Numerical  integration  of  E  from 
equations  (8)  to  (13)  was  therefore  used  to  produce  the  curves 
of  Fig.  9.  They  snow  the  view  factor  F(y#3)  for  the  sphere  as  a 
function  of  sensor  angle  for  various  target  heights.  As  before, 
the  curves  are  replotted  in  Fig.  10  normalized  to  unity  at  Y=0r 
to  show  the  angular  depen  lence  more  clearly.  Figure  11  shows 
F(y,B)  as  a  function  of  height  for  several  values  of  y. 

Special  Cases  When  the  target  height  is  either  very  small  or  very 

large,  the  formulae  of  equations  (15)  and  (16)  can  be  used  to 
yield  cl  sed  form  expressions  for  the  view  factor  of  a  diffuse 
target,  since  they  are  easily  integrable.  The  resulting  expressions 
are 

F(y,6+|)  =  i  +  j  cos  y  (22) 
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SENSOR  ANGLE  (dtg) 


Fig.  9.  view  factor  of  diffuse  sphere  as  a  function 
of  sensor  angle  for  various  target  heights. 
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SENSOR  ANGLE  (deg) 


Fig.  10.  View  factor  of  diffuse  sphere  normalized  to  unity 
for  sensor  angle  Y=0;  various  target  heights. 


1.0 


Fig.  11.  view  factor  of  diffuse  sphere  as  a  function  of 
target  height  for  various  sensor  angles. 


(23) 


F<Y,j3-H))  » 
These  curves  are  shown  as 


_2 

3ir 


h=0 


sin  y  +  ( 7t~y )  cos 
and  h=°°  on  Figs. 


Yj  sin2 6 
9  and  10. 
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SUMMARY 


Simple  expressions  have  been  obtained  for  the  distribution 
of  earthshine  incident  on  a  spherical  target,  as  a  function  of 
target  height.  These  results  apply  to  the  case  of  a  uniformly 
and  isotropically  emitting  earth. 

Curves  for  the  reflected  signature  (or  view  factor)  measured 
by  a  sensor  at  an  arbitrary  location  have  been  obtained  by  inte¬ 
grating  the  reflected  earthshine  over  the  portion  of  the  target 
visible  to  the  sensor,  for  the  case  of  a  diffuse  target.  A  simple 
expression  is  given  for  the  view  factor  of  a  specular  target. 

The  view  factor  results  are  also  valid  for  clouds  of  small 

*» 

randomly  shaped  convex  particles  of  corresponding  surface 
characteristics . 
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